Dynamical critical phenomena in driven-dissipative systems 
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We explore the nature of the Bose condensation transition in driven open quantum systems, such 
as exciton-polariton condensates. Using a functional renormalization group approach formulated in 
the Keldysh framework, we characterize the dynamical critical behavior that governs decoherence 
and an effective thermalization of the low frequency dynamics. We identify a critical exponent 
special to the driven system, showing that it defines a new dynamical universality class. Hence 
critical points in driven systems lie beyond the standard dynamical classification of equilibrium 
phase transitions. We show how the new critical exponent can be probed in experiments with 
driven cold atomic systems and exciton-polariton condensates. 



Recent years have seen major advances in the explo- 
ration of many-body systems in which matter is strongly 
coupled to light [1] . Such systems include for example po- 
lariton condensates [2] , superconducting circuits coupled 
to microwave resonators [3, 4], cavity quantum electro- 
dynamics [5] as well as ultracold atoms coupled to high 
finesse optical cavities [6]. As in traditional quantum 
optics settings, these experiments are subject to losses, 
which may be compensated by continuous drive, yet they 
retain the many-body character of condensed matter. 
This combination of ingredients from atomic physics and 
quantum optics in a many-body context defines a qual- 
itatively new class of quantum matter far from thermal 
equilibrium. An intriguing question from the theoretical 
perspective is what new universal behavior can emerge 
under such conditions. 

A case in point are exciton-polariton condensates. Po- 
laritons are short lived optical excitations in semicon- 
ductor quantum wells. Continuous pumping is required 
to maintain their population in steady state. But in 
spite of the non-equilibrium conditions, experiments have 
demonstrated Bose condensation [2] and, more recently, 
have even observed the establishment of a critical phase 
with power-law correlations in a two dimensional system 
below a presumed Kosterlitz-Thouless phase transition 
[7]. At a fundamental level however there is no under- 
standing of the condensation transition in the presence of 
loss and external drive, and more generally of continuous 
phase transitions under such conditions. 

In this letter we develop a theory of dynamical critical 
phenomena in driven-dissipative systems in three dimen- 
sions. Motivated by the experiments described above we 
focus on the case of Bose condensation. It is a com- 
mon belief that such systems develop an effective tem- 
perature in steady state, and that phase transitions in 
driven-dissipative systems mimic in all aspects classical 
dynamical critical behavior at finite temperature. In cer- 
tain cases it is observed that an effective temperature in- 



deed develops [8-11]. We show below that while static 
properties of the driven-critical condensate are indeed 
governed by an effective temperature, the dynamical be- 
havior belongs to a different universality class. To this 
end, we relate and highlight the differences between the 
dynamics of a driven system near a critical point, and 
the well-known classification of dynamical critical phe- 
nomena around equilibrium phase transitions [12]. In 
particular, we identify a measurable universal exponent 
that betrays the driven nature of the system and would 
not be obtained in the dynamic response at equilibrium 
critical points. 

Open system dynamics- A microscopic description of 
driven open systems typically starts from a Markovian 
quantum master equation or an equivalent Keldysh ac- 
tion. However, the novel aspects in the critical dynamics 
of driven dissipative systems discussed below can be most 
simply illustrated by considering an effective phenomeno- 
logical description of the order parameter dynamics using 
a stochastic Gross-Pitaevskii equation 
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As we show below, this equation can be rigorously de- 
rived from a fully quantum microscopic description of 
the condensate when including only the relevant terms 
near the critical point. The different terms in (1) have 
a clear physical origin. X = lp ~ li is the effective gain, 
which combines the incoherent pump field minus the local 
single-particle loss terms, k, X > are respectively two- 
body loss and interaction parameters. The diffusion term 
D is not contained in the original microscopic model, but 
is generated upon integrating out high frequency modes. 
We therefore include it at the outset with a phenomeno- 
logical coefficient. Finally £ is a Gaussian white noise 
with correlations (C*(t,x)C(t',x')) = jS(t - f)S(x - x'). 
Such noise is necessarily induced by the losses, and sud- 
den appearances of particles due to pumping. 
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Without the noise term, Eq. (1) is similar to the dissi- 
pative Gross-Pitaevskii (dGP) equation previously used 
as a mean field description of the exciton-polariton con- 
densates [13, 14]. The dGP admits a stationary conden- 
sate solution with density \ip\ 2 = — x/k for x < 0, under- 
going a mean field transition to the vaccum state when 
X crosses zero. The chemical potential of the condensate 
/i = A|?/»| 2 is determined by dynamical stability [15]. We 
also note that the dGP equation used in [13, 14] does not 
feature the diffusion term D mentioned above, which is 
expected to play a role near the critical point. 

The noise term in Eq. (1) acts similarly to a temper- 
ature, in that it can drive a transition at finite particle 
density, thereby inducing critical fluctuations. As the 
equation of motion is cast in Langevin form, one might 
suspect that it can be categorized into one of the well- 
known models of dynamical critical phenomena classified 
by Hohenberg and Halperin [12]. However, this is not 
true in general. The models in Ref. [12] are constrained 
to have a very specific relation between the reversible and 
dissipative terms to ensure a thermal Gibbs distribution 
of the fields in steady state [16, 17]. This is equivalent 
to the requirement that both the coherent and dissipa- 
tive evolutions are generated by the same Hamiltonian. 
In the driven open systems we consider, there is no such 
requirement. The dissipative dynamics is determined by 
the intensity of the pump and loss terms, independently 
of the intrinsic Hamiltonian dynamics of the system. For 
the critical dynamics this implies that the scaling of the 
dissipative coupling constants is not necessarily tied to 
that of the reversible ones by the equilibrium constraint, 
opening up the possibility for novel dynamical universal- 
ity classes. 

Microscopic Model - Having illustrated the nature of 
the problem with the effective classical equation (1) we 
turn to a fully quantum description within the Keldysh 
framework. Our starting point is the Markovian dissipa- 
tive Keldysh action [18] 
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S q are the so-called classical and quantum fields, 
given by the symmetric and anti-symmetric combinations 
of the fields on the forward and backward parts of the 
Keldysh contour [19, 20]. The inverse Green's functions 
are given by P B = id t + AW 2 + (i + i\, P A = P R \ 
P K = ij. 

Vanishing of the mass scale x defines a Gaussian fixed 
point with dynamical critical exponent z — 2 (w ~ k z ). 
Canonical power counting determines the scaling dimen- 
sions of the fields and interaction constants with respect 
to this fixed point: at criticality, the spectral components 
of the Gaussian action scale as P R / A ~ k 2 , while the 



Keldysh component generically takes a constant value, 
i.e., P K ~ k°. Hence, to maintain scale invariance of 
the quadratic action, the scaling dimensions of the fields 
must be [<p c ] = ^ and [<fi q ] = ^±2. From this result 
we read off the canonical scaling dimensions of the in- 
teraction constants; in particular, for a local interaction 
vertex containing n c classical and n q quantum fields, we 
find [An cl n ] = d + 2 — n c ^=^ — riq^ 2 -. Hence, the upper 
critical dimension for the quartic couplings is d = 4 as 
in finite temperature thermodynamic equilibrium. More- 
over, in the case of interest d = 3, local vertices contain- 
ing more than two quantum fields (n q > 1 is required by 
causality [19, 20]) or more than five classical fields are 
irrelevant. The only marginal term with n q — 2 is the 
Keldysh component of the single-particle inverse Green's 
function, i.e., the noise vertex. In this sense, the criti- 
cal theory is equivalent to a stochastic classical problem 
[21, 22], as previously observed in [8, 23]. But as noted 
above it cannot be a priori categorized in one of the dy- 
namical universality classes [12] subject to an intrinsic 
equilibrium constraint. 

Functional RG - To focus on the critical behavior we 
turn to a functional RG formulation for the field the- 
ory (2). The original equilibrium formulation due to Wet- 
terich [24] was adapted to the Keldysh real time frame- 
work in [25, 26]. Central to this approach is the func- 
tional T\[(f) c ,(p q ] defined by [27] 



(3) 

Here A6>a is a regulator function which suppresses con- 
tributions to the above path integral from modes with 
spatial wave-vector below the running cutoff A (see SI 
for details). Thus Ta interpolates between the classical 
action S, when A equals the UV cutoff A , and the ef- 
fective action functional T[(j> c , <j> q ] [28] when A — > 0. The 
latter includes the effects of fluctuations on all scales. 
The equation 
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describes the flow of the interpolating functional as a 

(2) 

function of the running cutoff A. Here T\ and R A de- 
note the second field variations of Ta and AS\ respec- 
tively, which makes Eq. (4) a functional differential equa- 
tion. For the description of general equilibrium [29] and 
Ising dynamical [30] critical behavior the functional RG 
gave results that are competitive with high-order epsilon 
expansion and with Monte Carlo simulations already in 
rather simple approximation schemes. 

We solve Eq. (4) approximately using an ansatz for T\ 
which includes all couplings that are relevant or marginal 
at tree level. Specifically, we take 
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The fact that the spectral components of the effective ac- 
tion depend only linearly on <p q allowed us to introduce an 
effective potential U determined by the static couplings. 
U (p c ) = \u (p c — po) 2 + g-O' (p c — po) 3 is a function of the 
t/(l) invariant combination of classical fields p c — <j)%<j)c 
alone. It has a mexican hat structure ensuring dynami- 
cal stability. With this choice we approach the transition 
from the ordered side, taking the limit of the stationary 
state condensate po = 4>l4> c \ s t = <f>Q<pQ —t 0. 

All the parameters appearing in (5) including the sta- 
tionary condensate density po are functions of the run- 
ning cutoff A. Hence, the functional flow of Ta is reduced 
by means of the approximate ansatz to the flow of a finite 
number of coupling constants g = (Z, K, po, u, v! , 7) 
determined by the /3-functions Ac*Ag = /3 g (g) (see SI). 
The critical system is described by a scaling solution to 
these flow equations. In practice, it is obtained as a fixed 
point of the flow of dimensionless (according to the above 
power counting) renormalized couplings, which we derive 
in the following. First we rescale couplings with Z , 

K = K/Z, u = u/Z, v! = u'/Z, 7 = 7/ \Z\ 2 . (6) 

Coherent and dissipative processes are encoded, respec- 
tively, in the real and imaginary parts of the renormal- 
ized kinetic coefficient K = A + iD and the two- and 
three-body couplings u = A + in and u' = A' + in' . The 
relative strength of these contributions is thus measured 
by the (dimensionless) ratios tk — A/D, r u — X/k, and 
?V = A'/k', and we will find it convenient to work with 
these ratios instead of the respective real parts. The two- 
and three-body loss coefficients k and n' and the conden- 
sate density po become dimensionless by means of the 
rescalings 
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The flow equations for the couplings r = (tk, t u , r u >) 
and s = (w, k, k') T form a closed set, 

A9 A r = &(r, s), A0 A s = /3 s (r, s) (8) 

(see SI for the explicit form). As a consequence of the 
transformations (6) and (7), these /3-functions acquire 
a contribution from the running anomalous dimensions 
rj a (r, s) = — A<9a In a associated with a = Z,D, 7. 

Critical properties - The universal behavior near the 
critical point is controlled by the infrared flow to a 
Wilson-Fisher like fixed point. The values of the cou- 
pling constants at the fixed point, determined by solving 
/3 s (r*,s*) = and /3 r (r», s*) = 0, are given by: 

I"* = {rx* 1 Tu* 5 Til'* ) = 0, 



(w*,k*,k*) ~ (0.475,5.308,51.383). 



(9) 



The first line implies that the fixed point action is purely 
imaginary (or dissipative), as in Model A of Hohenberg 
and Halperin [12], cf. Fig. 1 (c). We interpret the fact 
that the ratios of coherent vs. dissipative couplings are 
zero at the fixed point as a manifestation of decoherence 
at low frequencies in an RG framework. The coupling 
values are identical to those obtained in an equilibrium 
classical 0(2) model from functional RG calculations at 
the same level of truncation [29] . 

Let us turn to the linearized flow, which determines 
the universal behavior in the vicinity of the fixed point. 
We find that the two sectors corresponding to s and r 
decouple in this regime, giving rise to a block diagonal 
stability matrix 



d 
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(10) 



where Sr = r, 5s = s — s*, and N, S are 3x3 matrices 
(see SI). 

The anomalous dimensions entering this flow are found 
by plugging the fixed point values r*, into the expres- 
sions for rj a (r, s). We obtain the scaling relation between 
the anomalous dimensions r\z = — ^7, valid in the uni- 
versal infrared regime. This leads to cancellation of r\z 
with T] 7 in the static sector S (see SI) . The critical prop- 
erties in this sector, encoded in the eigenvalues of S, be- 
come identical to those of the standard 0(2) transition. 
This includes the correlation length exponent v « 0.716 
and the anomalous dimension 77 « 0.039 associated with 
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initial flow 



Re 

linearized IR flow 



Re 



fixed point 



FIG. 1. Flow in the complex plane of dimensionless renor- 
malized couplings, (a) The microscopic action determines the 
initial values of the flow. Typically, the coherent propagation 
will dominate over the diffusion, A D, while two-body col- 
lisions and two-body loss are on the same order of magnitude, 
A ~ k, with a similar relation for the marginal complex cou- 
pling u' . The initial flow is non-universal, (b) At criticality, 
the infrared (IR) flow approaches a universal (independent 
of the initial conditions) linear domain encoding the critical 
exponents and anomalous dimensions, (c) The Wilson-Fisher 
fixed point describing the interacting critical system is purely 
imaginary. 
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the bare kinetic coefficient K. These vaiues are in good 
agreement with more sophisticated approximations [31]. 

The equilibrium-like behavior in the S sector can be 
seen as a result of an emergent symmetry. Locking of 
the noise to the dynamical term implied by the scal- 
ing relation rjz = — ?7 7 leads to invariance of the long 
wavelength effective action (times i) under the trans- 
formation $ c (t,x) -> $ c (-t,x),<I> g (i,x) -S- <f> g (-t,x) + 
|o-*3t$ c (-t,x),i -> -i with $„ = (<j> v , 0* ) T , v = (c,q), 
o~ z the Pauli matrix. This symmetry is a direct mani- 
festation of a classical fluctuation-dissipation relation. It 
generalizes the symmetry noted in Refs. [32, 33] to mod- 
els that include also reversible couplings compatible with 
an equilibrium stationary state. Unlike such Hohcnbcrg- 
Halperin type models here the symmetry is not imposed 
at some mesoscopic level of the theory, but rather is emer- 
gent at the critical point. We interpret this finding as an 
effective low-frequency thcrmalization mechanism of the 
driven system at criticality. 

The key new element in the driven-dissipative dynam- 
ics is encoded in the decoupled "drive" sector (the 3x3 
matrix ./V in our case). It describes the flow towards 
the emergent purely dissipative Model A fixed point (see 
Fig. 1 (b)) and thus reflects a mechanism of low frequency 
decoherence. This sector has no counterpart in the stan- 
dard framework of dynamical critical phenomena and is 
special to driven-dissipative systems. In the deep infrared 
regime, only the lowest eigenvalue of this matrix governs 
the flow of the ratios. This means that only one new 
critical exponent r\ r rj —0.101 is encoded in this sector; 
below we indicate how it affects single particle observ- 
ables. Just as the dynamical critical exponent z is inde- 
pendent of the static ones, the block diagonal structure 
of the stability matrix ensures that the drive exponent is 
independent of the exponents of the other sectors. 

The fact that the inverse Green's function in Eq. (5) is 
specified by three real parameters, KeK,lmK, and \Z\ 
(the phase of Z can be absorbed by a U(l) transforma- 
tion) allows for only the three independent anomalous 
dimensions: r/o, Vz and the new exponent rj r . Hence 
the extension of critical dynamics described here is max- 
imal, i. e., no further independent exponent will be found. 
Moreover this extension of the purely relaxational (Model 
A) dynamics leads to different universality than an ex- 
tension that adds reversible couplings compatible with 
relaxation towards a Gibbs ensemble. The latter is ob- 
tained by adding real couplings to the imaginary ones 
with the same ratio of real to imaginary parts for all 
couplings [34-37]. Such an extension adds only an inde- 
pendent lxl sector N to the purely relaxational problem, 
for which we find r/R = —0.143 ^ r) r . This proves that 
the independence of dissipative and coherent dynamics 
defines indeed a new non-equilibrium universality class 
with no equilibrium counterpart. 

Experimental detection - The novel anomalous dimen- 
sion identified here leaves a clear fingerprint in single- 



particle observables accessible with current experimen- 
tal technologies on different platforms. Using the RG 
scaling behavior of the diffusion and propagation coef- 
ficients D ~ hr 7 ^, A = Dr K ~ h-Vr-m^ we ob- 
tain the anomalous scaling of the frequency and mo- 
mentum resolved, renormalizcd retarded Green's func- 
tion G*(w,q) = (w - A |q| 2 -'™ + iD |qr" D )-\ 
with Aq and Dq non-universal constants. The peak 
position and width implied by the complex dispersion 
u! Ri Ao |q| 2 ' 22 — iDo |q| 2 12 can be measured in ultra- 
cold atomic systems by RF spectroscopy [38] . In exciton- 
polariton condensates, the dispersion relation can be re- 
constructed using homodyne detection techniques [39], 
which allows one to resolve both real and imaginary part 
of the retarded response separately. 

Finally, we estimate the extent of the universal critical 
domain governed by the linearized regime of the Wilson- 
Fisher fixed point, which provides us with an estimate 
of the energy resolution \c necessary to probe the crit- 
ical behavior. This is done by calculating the Ginzburg 
scale, i.e., the distance from the phase transition where 
fluctuations become dominant [28]: We equate the bare 
distance from the phase transition \ to the correspond- 
ing one-loop correction, yielding \c = 8njn/ D 3 / 2 . The 
diffusion constant D can be obtained from experiment 
by the same techniques as already indicated. A two-loop 
estimate gives the scaling D ~ A 2 n 2 , K 2 n 2 . 

Conclusions - We have developed a Keldysh field the- 
oretical approach to characterize the critical behavior of 
driven-dissipative three dimensional Bose condensates at 
the condensation transition. The main result presents 
a hierarchical extension of classical critical phenomena. 
First, all static aspects are identical to the classical 
0(2) critical point. In the next shell of the hierarchy 
a sub-class of the dynamical phenomena is identical to 
the purely dissipative model A dynamics of the equilib- 
rium critical point. Finally we identify manifestly non- 
equilibrium features of the critical dynamics, encoded 
in a new independent critical exponent that betrays the 
driven nature of the system. 

It is an interesting question for future work how the 
above conclusions are modified when applied to systems 
of lower dimensionality. Particularly pertinent for un- 
derstanding current experiments with exciton-polariton 
condensates in semiconductor quantum wells is the case 
of two dimensional systems. 
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FUNCTIONAL RENORMALIZATION GROUP 
EQUATION 

Our approach to studying critical dynamics is based on 
the Wetterich functional renormalization group equation 
[1] adapted to the Keldysh framework 



9 A r A = - Tr 



(rf + r a ) 'oaRa 



(SI-1) 



In the first part of this supplement we discuss the objects 
appearing in (SI-1), namely the second functional deriva- 

(2) 

tive and the cutoff function Ra- Next we explain how 
a closed set of flow equations for a finite number of cou- 
pling constants can be obtained from the functional flow 
equation. Finally we detail the linearized equations for 
the infrared flow to the Wilson-Fisher fixed point from 
which the critical properties are inferred. 

In suitable truncation schemes, results from high order 
epsilon expansion can be reproduced from the exact flow 
equation (SI-1). In our practical calculation, we approach 
the critical point from the ordered phase. This allows us 
to calculate the anomalous dimensions at one-loop order, 
due to the presence of a finite condensate during the flow. 
Results obtained in this way have proven to be compet- 
itive with high-order epsilon expansion or Monte Carlo 
simulations, as referenced in the main text. 



THE SECOND VARIATIONAL DERIVATIVE 

(2) 

The second variation T A with respect to the fields is 
the full inverse Green's function at the scale A, which 
in the case of an interacting theory is field dependent. 
Practically we work in a basis of real fields, related to 
the complex fields by 



x*AQ) 



i/ii 



MQ) 



(SI-2) 



where v = c, q is the Keldysh index. We gather the result- 
ing four independent field components in a field vector, 

x(Q) = (xc,i(Q),XcAQ),x q AQ),x q AQ)f ■ (si-3) 



In this basis, T^ 2 ' is defined as 



(rf). (q,q') 



6 2 T* 



(SI-4) 



a 8 Xi (-Q)SxAQ') 
which is a matrix in the discrete field index i = 1, 2, 3, 4 
and in the continuous momentum variable Q = (cj, q) col- 
lecting frequency and spatial momentum. Accordingly, 
the trace in (SI-1) involves both an integration over mo- 
menta and a sum over internal indices. 

(2) 

T A (Q, Q') is conveniently decomposed into a constant 
part and a fluctuation part. The latter is a polynomial 
in momentum-dependent fields and, therefore, a non- 
diagonal matrix in momentum space. In contrast, the 
constant part is obtained by (i) inserting spatially con- 
stant field configurations, i.e., x(Q) = X${Q) m momen- 
tum space, and (ii) evaluating them at their equilibrium 
values in the ordered phase. These read 

x(Q)| cq = (V2po, 0,0, 0) T <5(Q). (SI-5) 

(Without loss of generality we choose the condensate am- 
plitude to be real.) As a result, the constant part is di- 
agonal in momentum space, 



Pa{Q)5{Q-Q') = V { 1\q,Q' 



I cq 7 



(SI-6) 



and is structured into retarded, advanced, and Keldysh 
blocks, 



PAQ) 



p a (Q) 

pK 



(SI-7) 





(For notational simplicity, we suppress the scale index 
A for the different blocks and their respective entries.) 
The retarded and advanced blocks are mutually hermi- 
tian conjugate (we decompose Z and K into real and 
imaginary parts, Z = Zr + iZj, K = A + iD), 

-iZiu - Aq 2 - 2Rc(u)po iZ R u> - Dq 2 
- 2 Im(u)p -iZjOJ - Aq 2 



P iQ) - 1 -iZ R u; + Dq 2 

P A {Q) = (P K (Q)) t . 
Note that dct P R (Q = 0) 



(SI-8) 
0; the 



det P A (Q = 0) 
existence of a gapless mode associated to the broken U(l) 
symmetry is thus ensured in our truncation at all scales 
A. For the Keldysh block we have 

P K = *7l. (SI-9) 
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THE REGULATOR FUNCTION 

The cutoff ASa is used in Eq. (3) of the main text to 
generate the effective action Ta from the microscopic ac- 
tion S by suppressing contributions from momenta below 

(2) 

A. Its second functional derivative Ra — AS A enters 
the exact flow equation (SI-1). We choose an optimized 
cutoff function [6] of the form 



Ra(Q) = (q 2 - A 2 ) 9(A 2 - q 2 ) 



R R 
R A 



(SI-10) 



where 



R H = 



(SI-ll) 



Due to the 6>-function in (SI-10), in the regularized in- 
verse Green's function 



Pa + Ra, 



momenta q 2 smaller than the running scale A 2 acquire 
an effective mass oc A 2 and we have det G^ 1 (Q = 0) ^ 0, 
which ensures that momentum integrals over Green's 
functions arc infrared convergent. Note that it is suf- 
ficient for Ra to modify only the retarded and advanced 
blocks (i.e., the spectrum) of the inverse Green's func- 
tion. The choice of a frequency-independent cutoff allows 
us to perform frequency integrals analytically. 

The interpolation property of Y a between the classical 
action S and the effective action T is guaranteed by the 
limiting behavior [7] 



lim Ra 



A; 



o- 



lim Ra 

A 2 ->0 



0. 



(SI-13) 



FLOW OF THE EFFECTIVE POTENTIAL 

In equilibrium problems, an important object for prac- 
tical calculations is the effective potential. It describes 
the homogeneous part of the effective action and is ob- 
tained by evaluating the full effective action at spatially 
homogeneous field configurations, U — ^/^\ x ^q)— x s(q) 
(fl is the quantization volume). In the framework of a 
derivative expansion, a closed flow equation can be de- 
rived for this object, which serves as a compact generat- 
ing functional for the flow of all local couplings to arbi- 
trarily high order. Here we provide the Keldysh analog 
of this construction, where the key difference roots in the 
occurrence of two field variables (f> c ,<f> q , in contrast to a 
single field in equilibrium. However, for a theory which 
obeys the power counting discussed in the main text, we 
can parameterize the homogeneous part of the effective 
action as 



V 



dU 

d<j) c ' 



+ (SI-14) 



with U = U(4>*<p c ) dependent on the U(l) invariant com- 
bination of classical fields only, this function thus be- 
ing the direct counterpart of the effective potential. A 
flow equation can be derived for the auxiliary object V, 
which reads (we introduce a dimensionless scale deriva- 
tive de = Ad a) 



d e V 



Q 



tr [g A (Q)d t R\{Q)]. 



(SI-15) 



Here, the inverse of Qa is obtained from the full sec- 
ond functional variation by evaluating it at homogeneous 
field configurations (step (i) above Eq. (SI-5)), how- 
ever without inserting the equilibrium values (step (ii)): 
G^ 1 = t a 2) \ x (q) =x s(q) + Ra- Qa is then diagonal in mo- 
mentum space, and so the trace over momentum space in 
Eq. (SI-1) reduces to a single momentum integration, giv- 
ing rise to the above compact form. In contrast to G^ 1 , 
(SI-12) ^A nas a non-vanishing upper left block P H . However, 
it vanishes when the background fields are set to their 
equilibrium values, P H | = 0, which is a manifestation 
of causality in the Keldysh formalism [4, 5]. 



From this equation we obtain the /3-functions for the 
momentum-independent couplings by evaluating appro- 
priate derivatives with respect to the U(l) invariants 



Pcq 



J qci 



Pq 



V (SI-16) 



at their equilibrium values p c | oq = p , p cq \ cq = p qc \cq = 
Pq\cq — 0. Specifically, we use the projection prescrip- 
tions 



1 



depo 


= /3 P0 


den 


= Pu 


d f u' 


= Pu> 


da 


= /?7 



(SI-17) 



l J cql->qc 



d,V 



cq 



Calculation of the explicit expressions here and below is 
largely automatized using Mathematica. 



FLOW OF THE INVERSE PROPAGATOR 

While the flow equation for the effective poten- 
tial (SI-15) generates /3-functions for all momentum- 
independent couplings, the flow of the complex dy- 
namic Z and kinetic K couplings, which constitute the 
momentum-dependent part of the effective action (5), is 
determined by the flow equation for the inverse propa- 
gator. We obtain the latter by taking the second vari- 
ational derivative of the exact flow equation (SI-1) and 
setting the background fields to their equilibrium values 
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Eq. (SI-5), 

dtPA,ij(Q) = 

% -(tx [G 2 A (Q' - Q)diR\(Q' - 0)7,G A (Q) 7 , 

1 JQ> 



+G A (Q' - 0)7,G|(Q)9,i?A(Q) 7j ] , (SI-18) 



where 



liS 5(P -P' + Q) 



5T 



SxiiQ) 



(SI-19) 



cq 



In Eq. (SI-18) we omit tadpole contributions oc T*- 4 ' 



A > 



which do not depend on the external momentum Q and 
hence do not contribute to the flow of Z or K. For these 
we use the projection prescriptions 



8 t Z = fi z = tr [(1 + tr y ) d t P R {Q)] 

d e K = fi R = a q2 [d e P 2 R 2 (Q) + id t P&{Qj\ 



Q=0 



Q=0 



(SI-20) 



The /?-functions (SI- 17) and (SI-20) constitute the com- 
ponents of /? g = (Pz,(3K,Pp ,Pu,i3u>,fa) T - 



RESCALED FLOW EQUATIONS 

We write the flow equation for the complex dynamic 
coupling Z in the form 



a e z = -vzZ. 



(SI-21) 



The anomalous dimension rj z is an algebraic function of 
the rescaled couplings (6) and po- The same applies to 
the /3-functions of the latter, 

1 



d e K = (3 K = rj z K+-P R , 



diu 



d t u' = fi u , = rjzu' + gPv, 
dil = P-y = (vz + V*z)l+ T7^i$i- 



(SI-22) 



In particular, the very right expressions in these equa- 
tions (Pk/Z etc.) are functions of the rescaled couplings 
alone. In terms of these variables, therefore, all explicit 
reference to the running coupling Z is gone, and we have 
effectively traded the differential flow equation for Z for 
the algebraic expression for its anomalous dimension rj Z - 
All couplings except for 7 are complex valued. Taking 
real and imaginary parts of the /3-functions for K, u, and 
u' yields the flow equations for A, D, A, k, A', and k' 
respectively, 

d t A — Pa — Re Pk 7 d e D = /3 D =Im/3 K , 
8 e \ = p x = Re(3 u , d e n = (3 K = lm(3 u , (SI-23) 
d e X' = fi\> = Re/3 U >, den' = fi K < = Im/3„/. 



The /^-functions for the ratios r = (r K , r u , r u > ) are then 



d f r K = f3 rK = —f3 A - /?£>, 



dtr u = f3 ru = —f3\ - —p K , 

K K 

der u , = (3 ru , = — fix' - ^fi*'- 



(SI-24) 



(SI-25) 



K K 

The number of flow equations can be further reduced by 
introducing anomalous dimensions for D and 7, 

deD = -r] D D, 

da = -va- 

As for the dynamic coupling Z in terms of the rescaled 
variables K, u, vl , 7 and po, all explicit reference to D 
and 7 drops out, and we obtain for the couplings s = 
(w, Y defined in Eq. (7) 

d e w = fi w = - (2 - 7] D ) w + ^/3 K + ^t^po: 

8ek = fi k = -(l-2r lD +r ]l )k+^j^fi K , (SI-26) 



deh 1 = fix = - {-'irj D + 2?7 7 ) k' 



4D 



;fi K >. 



In summary, the transformations (6) and (7) result in the 
closed system (8) for r and s with fi r = (fi rK , (3 ru , fi Tu , ) 
given by Eq. (SI-24) and fi s = (fi w , fi k , fi k ,) T given by 
Eq. (SI-26). The flows of Z, D, and 7 are decoupled 
and determined by the anomalous dimensions (SI-21) 
and (SI-25), which arc themselves functions of r and s. 



CRITICAL PROPERTIES 

For the analysis of critical behavior, we need to find a 
scaling solution to the flow equations for the bare cou- 
plings or, equivalently, a fixed point r*, of the flow of 
dimensionless rescaled couplings, 



/? r (r*,s*) = /3 s (r*,s*) = 0. 



(SI-27) 



This non-linear algebraic set of equations has a non- 
trivial solution given by Eq. (9). In order to characterize 
the infrared flow in the vicinity of the fixed point (en- 
coding the critical exponents we are interested in here), 
we study the flow of the couplings linearized around the 
fixed point, cf. Eq. (10). The stability matrices N and S 
in this equation read explicitly 



= Vjfir 



s = v!fi s 



0.0525 
-0.0002 
0.4976 

-1.6204 
-3.1828 
-15.3743 



0.0586 0.0317\ 
-0.0526 0.1956 , 
-2.3273 1.9725/ 

(SI-28) 

0.0881 0.0046\ 
0.2899 0.0363 , 
-42.2487 2.1828/ 
(SI-29) 
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without coupling between r and s sectors. At present 
we cannot rule out that an extended truncation would 
couple them. However, since we already include all rel- 
evant and marginal couplings, we expect the decoupling 
to be robust or at least approximately valid to a good 
accuracy. 

The infrared flow of Z, D, and 7 is determined by 
the values of the respective anomalous dimensions at the 
fixed point. Equations (SI-21) and (SI-25) imply the scal- 
ing behavior 

Z~A-" Z , L>~A-" D , 7~A _ ^ (SI-30) 

for A — > 0. While t)d and ?y 7 describe the flow of real 
quantities and are, therefore, themselves real by defini- 
tion, r\z is in general a complex valued function of r and 
s. At the fixed point, however, the imaginary part van- 
ishes, 

lmjj z = 0, (SI-31) 

which ensures scale invariance of the full effective action 
at the critical point. 

As is indicated in the main text, the emergence of 0(2) 
model critical properties in the sector s is due to the scal- 
ing relation r\z = —ri-y, which ensures that these anoma- 
lous dimensions compensate each other in the /3-functions 
for the couplings s. This can be seen most simply by ex- 



pressing, e.g., k in terms of bare quantities, 

= jhnju) 7 Im(u/Z) 
' 2AIm(A^) 2 2Alm(K/Zy [ ' ' 

In this form it is apparent that the scaling <~ A _r?z which 
applies to both Z and I/7 drops out. Similar arguments 
hold for w and k'. Alternatively, the cancellation of rjz 
and ?7 7 in the /3-functions can be seen explicitly by insert- 
ing Eqs. (SI-22) and (SI-23) in (SI-26). What remains is 
a dependence on 77 = r/D+i]z which is just the anomalous 
dimension associated with the bare kinetic coefficient K. 
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